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MATRIX ELEMENTS OE LIGHT QUARKS 

G. Martinelli’^ 


Dip. di Fisica, Univ. “La Sapienza” and INFN, Sezione di Roma, P.le A. Moro, 1-00185 Rome, Italy. 

Lattice results for matrix elements of light quark operators are reviewed. The discussion is focused on recent theoretical 
progress and new numerical calculations which appeared after the Lattice 2000 Conference. 


1. Introduction 

In this talk a review of theoretical and numerical 
results for matrix elements of light quark operators 
is presented. The review is focused on the defini¬ 
tion, renormalization and improvement of operators 
relevant to different processes and on the evaluation 
of the corresponding physical amplitudes. At this 
Conference the main topics covered in the parallel 
sessions were electromagnetic and weak form factors, 
structure functions, K^-K^ mixing, kaon decay am¬ 
plitudes and light-cone wave functions. Correspond¬ 
ingly the problem of the renormalization and im¬ 
provement of the operators was also discussed. Re¬ 
sults with perturbative and non-perturbative, gauge 
invariant or gauge dependent renormalization meth¬ 
ods, using Ward identities, with unimproved or 
improved operators, with Wilson like (improved, 
twisted, ...), Domain Wall or overlap fermions were 
presented. In total I counted 38 talks in the parallel 
sessions and 13 posters. It is clear that it would have 
been impossible to discuss all this material, to which 
two accurate and complete review talks by L. Lel- 
louch and S. Sint ||] were dedicated at Banga¬ 
lore. A selection of the subjects has been then un¬ 
avoidable and I apologize to the authors of submitted 
contributions and to the speakers of talks given in the 
parallel sessions which I had to omit in my review. I 
also warmly thank all collegues who have kindly sent 
material and information for the preparation of this 
talk. 

The most exciting news at this Conference are 
the intriguing results for e'/e presented by the CP- 
PACS ^ and RBC |Q| Collaborations. Thus I start 
from these results to discuss a few physics issues of 
general interest. The renormalization of lattice quark 
operators, the ultraviolet problem, and the calcu¬ 
lation of final state interaction effects, the infrared 
problem, which is essential for kaon decays, are the 
next subjects. This will give me the opportunity to 
show several results that were presented in paral¬ 
lel talks at this Conference. More physics problems 
which are interesting for the lattice comunity can be 
found in the talk by M. Beneke at this Conference S 
and by J. Donoghue at the Kaon 2001 Conference H. 


The main topics in this review are: i) the determi¬ 
nation of the AI =1/2 and AI = 3/2 AT ^ tttt am¬ 
plitudes ii) the main contributions to e'/e 

due to the electropenguin operators {Qr and Qs in 
the standard notation) and to the strong penguin op¬ 
erator Qq iii) matrix elements of SUSY op¬ 

erators which can be computed on the lattice ||] ; iv) 
studies of the chiral expansion of the relevant ampli¬ 
tudes in finite and infinite volumes, in the quenched 
and unquenched cases 0-11 ; v) studies of Final 
State Interactions (FSI) in lattice simulations i ; vi) 
the extraction of the physical amplitudes from ma¬ 
trix elements computed in the Euclidean space-time 
on a finite volume HU- The problem of the renor¬ 
malization and improvement of the relevant opera¬ 
tors will be mentioned and discussed along with the 
presentation of the different arguments. Although 
the discussion is mainly focused on kaon physics, 
many considerations are general and apply to the 
calculation of all matrix elements. 

2. General physics issues 

Flavour changing neutral currents (FCNC) have 
been and are the main source of information on quark 
masses, couplings and CP violation. It is sufficient to 
mention that they led to the GIM theory and that 
the large B-B mixing was the first strong indica¬ 
tion that the top quark was very heavy, long be¬ 
fore precision tests at LEP or its discovery. They 
are also a window on physics beyond the Standard 
Model (SM). Since in the SM FCNC cannot occur at 
tree level, their effects compete with contributions 
from other heavy sectors of the theory as the extra 
particles predicted, for example, by Supersymmetry. 
This explains the large experimental and theoretical 
efforts in FCNC studies. In table |^, a list of FCNC 
processes for which experimental measurements ex¬ 
ist is presented. In two cases the comparison of the 
theory with experiments requires the knowledge of 
the hadronic matrix elements. Lattice results for the 
matrix element of the AS = 2 operator relevant for e, 
encoded in the bag parameter Bk, exist since a long 
time. This year, for the first time, we have two com¬ 
plete calculations of the matrix elements of the oper¬ 
ators entering e'/e (and also ReAo and ReA 2 ) which 
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Table 1 

Experimental measurements and theoretical predictions for several neutral current processes: the CP violation parameters 
e and t jt, the neutral B meson mixing AMs and AMd, radiative inclusive & —> 57 decays and In the 

latter two cases the calculation can be done by using the HQET expansion and the experimental A" —> tt weak current 
matrix element respectively. Thus a lattice calculation of the hadronic amplitudes is not needed. 


Quantity 


Experiments 


Theory 


e 

e'/e 

AMs/AMd 

BR{B ^ Xsj) 
BR{K^ TT^vv) 

2.271 ± 0.017 X 10"^ iq 
17.2 ± 1.8 X lO-'^ |r^ 

> 30 (95%G.F.) 

3.11 ±0.39 X 10-'‘ m 

1 . 5 ;^^ X 10-1° @ 

q{l- p) Bk 
-7-l30 X 10-4 

3.50 ± 0.50 X 10 4 No Lattice Needed 

0.8 ± 0.3 X 10-4° No Lattice Needed 

Table 2 

Lattice results for AI = 1/2 transitions 
numbers are also given. 

using K TT matrix elements from RBC and CP-PACS. The experimental 

Reference 

FeAo 

H.e^2 

ReAa / ReA2 

e'/e 

CP-PACS @ 

Exps. [IT- 3 2001 

16-L 21 X 10-« 

29-L 31 X 10-® 
33.3 X 10-® 

1.3 -L 1.5 X 10"*^ 

1.1-L 1.2 X 10"® 
1.5 X 10"® 

9-l 12 
24-L27 

22.2 

-7-l -2 X 10-4 
- 8 -l -4 X 10-4 

17.2 ± 1.8 X 10-4 


can be compared with the experimental results. Here 
I want first to examine the physics implications of the 
RBC and CP-PACS calculations which are in total 
disagreement with the experimental measurements 
of e'/e by NA48 and KTeV, table ||. If confirmed, 
these results imply other contributions to e'/e from 
physics beyond the SM. Let me give an explicit ex¬ 
ample which could explain the discrepancy. In the 
effective Hamiltonian, besides the dimension-six four 
fermion operators considered in refs. there are 
two dimension-five chromomagnetic operators of the 
form 


ng = c; 
where 


Q 


C-Q- 


( 1 ) 


9 


g - ± -SRa^''edLG%) . ( 2 ) 

The chromomagnetic contribution exists in the SM 
(although it has not been evaluated by ||,^) but it is 
expected to give a rather small contribution to e'/e 
since it is chirally suppressed 0{m\)) |211. SUSY 


effects may enhance the chromomagnetic term | 22 [ | 
and this could be the explanation of the difference 
between lattice results and experiments. A large 
chromomagnetic term may trigger visible effects in 
hyperon decays and in CP violating K tttttt 
amplitudes Q . Its coupling is constrained also by e 
since it may give a large contribution there, although 
this is rather hard to estimate. Thus for the lat¬ 
tice comunity it is certainly very interesting to com¬ 
pute the relevant amplitudes of the chromomagnetic 
operators, similarly to what has been done for the 
electromagnetic one, which enters the CP violating 
Kr 7r°e+e“ decay p5|]. 


FCNC exist for different changes of flavour: 


(3) 


AF = 1, corresponding to e'/e, and AF = 2 pro¬ 
cesses, corresponding to e and B~B mixing, are the 
most studied. Although they are of fundamental 
importance for the strong CP problem, AF = 0 
FCNC processes have not received, instead, the nec¬ 
essary attention by the lattice community. AF = 0 
FCNC processes are also important because from 
them one may derive very stringent constraints for 
any extension of the SM. In this framework partic¬ 
ular important is the electric dipole moment of the 
neutron for which a very tight upper bound exists, 
djV' < 6.3 X 10“^® e cm [Q. After the first stud¬ 
ies at the end of the 80s | 2 ^, the calculation of the 
neutron electric dipole moment has been abandoned. 
This is very surprising since the renormalization of 
the relevant operators ot mqjsq) and the 

calculation of the necessary disconnected diagrams 
with stochastic sources is by now a common practice. 
There are several extra operators contributing to dj\f 
that appear in extension of the SM and the matrix 
elements of which have never been computed ||28| 


K 


AF=0 

eff 


= -'-C, qa^^^j^qFg, - '-Cc qG^ 


0^9 J abc^ pp^ V i 


(4) 


where Ce,c,g are the model dependent Wilson coef¬ 
ficients of the effective low energy theory. After this 
brief discussion of FCNC amplitudes which could 
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be computed on the lattice with interesting implica¬ 
tions, let us return to the results which have been 
presented at this Conference and to the progress 
which have been achieved since Lattice 2000. 

3. General framework 

Physical kaon weak decay amplitudes can be de¬ 
scribed with negligible error (of where 

fj, is the renormalization scale) in terms of ma¬ 
trix elements of the effective weak Hamiltonian 
(7r7r|7f^j^^ \K), written as combination of Wilson co¬ 
efficients and renormalized local operators 

■ (5) 

^ i 

The sum is over a complete set of operators, which 
depend on /i. In general there are 12 four-fermion 
operators and two dimension-five operators: a chro- 
momagnetic one and an electromagnetic one. As 
discussed before, in the SM the contribution of the 
dimension-five operators is usually neglected. The 
calculation of the matrix elements must be done non- 
perturbatively, and this is the role of the lattice, 
whereas the Wilson coefficients can be computed in 
perturbation theory. 

Wilson coefficients and renormalized opera¬ 
tors 

ForHf//=i the Wilson coefficients have been com¬ 
puted at the next-to-leading order |^]. The pertur¬ 
bative calculation is reliable provided that the scale 
^ is large enough, ^ ^qcd- In this respect cal¬ 
culations performed below the charm quark mass 
{rric ~ 1.3 GeV) are, in my opinion, suspicious. In 
fact, either ^ <C TOc, and then perturbation theory 
is questionable, or /r ~ me, and then the effective 
three-flavour weak Hamiltonian (with propagating 
up, down and strange quarks) cannot be properly 
matched to the four-flavour theory (up,down, strange 
and charm) because of the presence of operators 
of higher dimension which contribute at ©(/i^/m^). 
Wilson coefficients and matrix elements of the oper¬ 
ators Qi{n) separately depend on the choice of the 
renormalization scale and scheme. In the lattice ap¬ 
proach, it is possible to show that this dependence 
cancel in physical quantities, up to higher-order cor¬ 
rections in the perturbative expansion of the Wilson 
coefficients. 

Matching of bare (divergent) lattice operators, 
Qi{a) to the renormalized ones is obtained by com¬ 
puting suitable renormalization constants Zik (^a) 

"4*7=0,2(m) = (’’■7r|Q*(/7)|A:)7=o.2 

= ^ Zjfc(/7a)(7r7r|(5fc(a)|A:)/=o,2 , (6) 

k 


where a is the lattice spacing. The ultra-violet (UV) 
problem, which deals with the construction of fi¬ 
nite matrix elements of renormalized operators con¬ 
structed from the bare lattice ones, has been ad¬ 
dressed in a series of papers |3^-|3^ and is, at least in 
principle, completely solved. The remaining difficul¬ 
ties are practical ones. On the one hand, the A/ = 
1/2 operators suffer from power divergences in the 
ultraviolet cutoff, 1/a. These divergences, that can¬ 
not be subtracted using perturbation theory , can 
be eliminated by performing suitable numerical sub¬ 
tractions. The subtraction procedure, however, suf¬ 
fers from systematic uncertainties which are difficult 
to keep under control, see below. On the other hand, 
the perturbative calculation of the logarithmically di¬ 
vergent (or finite) Zik{^a) is rather inaccurate. Sev¬ 
eral non-perturbative methods have been developed 
in order to compute Ziki^ia) non-perturbatively [Q- 
-p^ and the uncertainties vary between 1%, for the 
simplest bilinear operators, to 10-;-25%, in the case of 
the four-fermion operators of interest. An accurate 
determination of Zn^^^a) for the A/ = 1/2 opera¬ 
tors has not been achieved to date and more work is 
needed in this direction. Let me also recall a method 
for a gauge-invariant non-perturbative renormaliza¬ 
tion of the relevant operators which was proposed 
in |4^ and never tried, in spite of the fact that it 
should be rather easy to implement. This proposal 
is based on the study, at short Euclidean distances 
\x\ <C I/Aqcd, of suitable correlators. For example, 
in order to renormalize bilinear operators of the form 
J{x) = q{x)Tq(x), one should compute numerically 
(0| J(a:)I J(0)|0) and compare the results with the ex¬ 
pression obtained in perturbation theory in a given 
popular renormalization scheme {MS or RI-MOM 
for example). The method avoids the problem of the 
window, a <C |a;| <C 1/Aqcd, by successive match¬ 
ing of the renormalization conditions and can 
be extended to four fermion operators. For bilinear 
operators, the necessary two-loop perturbative cal¬ 
culations, including the finite terms, have been com¬ 
pleted and the extension to four fermion operators is 
underway pl|| . 

Since numerical simulations are performed at fi¬ 
nite values of the lattice spacing, a~^ ^2-^4 GeV, 
another source of uncertainty in the determination of 
the matrix elements comes from discretization errors. 
They are of 0{Aqcdo), 0{\f^a) or 0{mqa), where 
p is a typical hadron momentum and mq the quark 
mass. The simplest strategy to reduce discretization 
effects consists in computing physical quantities at 
several values of the lattice spacing and then extrap¬ 
olate to a —> 0. A different approach is to reduce 
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discretization errors from 0{a) to O(a^) by improv¬ 
ing the lattice action and operators |39|] . A system¬ 
atic study of the improvement of four fermion oper¬ 
ators is still to be done. Note that with DWF 
or overlap Fermions the errors are automatically 
of O(a^) jiij. Discretization errors correspond to 
the matching problems in effective theories with a 
low cutoff recently discussed by Cirigliano, Donoghue 
and Golowich Q. In this respect the problem is 
softer for the lattice approach since numerical sim¬ 
ulations are already performed at relatively large 
scales. 

K ^ TT and K ^ tttt matrix elements 

Two main roads have been suggested in the past 
in order to obtain 2 (/^)- 

i) Compute the iF —> 0 and K ^ tt matrix ele¬ 
ments {0\Qi{ii)\K) and {TT\Qi{ii)\K) and then derive 
(7r7r|(5i(M)l^)/=o,2 using soft pion theorems @i3- 
In this case the K —>■ tttt amplitudes can be evalu¬ 
ated only at the lowest order of the chiral expansion 
and FSI effects in the physical amplitudes are lost. 
This problem poses serious difficulties to lattice cal¬ 
culations based on the extraction of the K —> tttt 
amplitudes from the K ^ tt matrix elements using 
chiral perturbation theory (yPT), as done by the CP- 
PACS 1^ and RBC Collaborations 0. The impor¬ 
tance of FSI effects for A/ =1/2 transitions was first 
noticed in refs. [jlTlJd^ ]. More recently it has been 
emphasized by Bertolini, Eeg and Fabbrichesi Hi 
and Pallante and Pich that they may have large 
effects also for e'/e. Although the quantitative eval¬ 
uation of FSI effects is still controversial there 
is a general consensus that, for a reliable calculation 
of kaon decay amplitudes, it is necessary to have a 
good theoretical control of FSI. 

ii) Compute directly (7r7r|Qj(/r)|A:)7^o.2 @0111 
1^ . The main difficulty in the latter case is due 
to the relation between K tttt matrix elements 
computed in a hnite Euclidean space-time volume 
and the corresponding physical amplitudes (the in¬ 
frared problem). The infrared problem arises from 
two sources: i) the unavoidable continuation of the 
theory to Euclidean space-time and ii) the use of 
a finite volume in numerical simulations. An im¬ 
portant progress has been achieved by Lellouch and 
Liischer (LL), who derived a relation between the 
lattice K —>■ tttt matrix elements in a finite volume 
and the physical kaon-decay amplitudes 0- An al¬ 
ternative discussion of boundary effects and the LL- 
formula, based on a study of the properties of corre¬ 
lators of local operators was presented in ref. . In 
this approach the LL-formula is derived for all elastic 
states under the inelastic threshold, with exponen¬ 


tial accuracy in the quantization volume. As a con¬ 
sequence the relation between finite-volume matrix 
elements and physical amplitudes, derived by LL for 
the lowest seven energy eigenstates, can be extended 
to all elastic states under the inelastic threshold. It 
can also been explicitly demonstrated how finite vol¬ 
ume correlators converge to the corresponding ones 
in infinite volume. The derivation of [p^ is based on 
the property of correlators of local operators which 
can be expressed, with exponential accuracy, both 
as a sum or as an integral over intermediate states, 
when considering volumes larger than the interaction 
radius and Euclidean times 0 < t ~ L. It can also 
been shown that it is possible to extract K tttt 
amplitudes also when the kaon mass, ttik, does not 
match the two-pion energy, namely when the inserted 
operator carries a non-zero energy-momentum. Such 
amplitudes are very useful, for example, in the deter¬ 
mination of the coefficients of operators appearing at 
higher orders in yPT, as illustrated by the numerical 
results for AI = 3/2 transitions presented in sec. 

Let us sketch now the derivation of the result with 
an illustrative example which is not explicitly written 
in 0] . Consider the following Euclidean T-products 
(correlators): 

G{t,tK) = {0\T[J{t)Q,i0)K\tK)]\0), 

G{t) = (0|T[J(t)J(0)]|0), 

GitK) = {0\T[KitK)K^0)\0), (7) 

where J is a scalar operator which excites (annhi- 

lates) zero angular momentum tttt states from (to) 
the vacuum and AT^ is a pseudoscalar source which 
excites a kaon from the vacuum {t > 0 and Ik < 0). 
At large time distances we have 

Git, Ik) —> 

V '^{Q\J\TnTin))v{'KTT{n)\Qi\K)v{K\K'^\Q)v 

n 

exp(-W„ t - niK tx), (8) 

Git) V y^(0| J|7r7r(n))y (7r7r(n)| J|0)v exp(—IF„ t) 

n 

GitK) ^ Vi{)\K\K)viK\K^\Q)veM-mKtK) 

Erom a study of the time dependence of G(t,<if), 
Gif) and Gitx) we may extract: 1) the kaon mass 
ruK', 2 ) the two-pion energies on the finite lattice 
volume, Wn- As shown by M. Liischer in |0|, the 
Wn are related to the infinite volume phase-shift of 
the two pions, 5(fc), via the following relations 

TT 
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n= 1,2,... , g= — , (9) 

where n is a non-negative integer and the func¬ 
tion is defined in |^. 3) the operator ma¬ 

trix elements on a finite volume {TnT{n)\Qi\K)v, 
{Tnr{n)\J\0)v and {K\K^0)v- Moreover by a suit¬ 
able choice of the lattice parameters it is possible to 
match one of the two-pion energies in such a way 
that rriK = ■ In practice, since it will be possi¬ 

ble to disentangle only the first few states, the matrix 
elements will be computed with n* = 0 -h 2 . 

The fundamental point is that it is possible to 
relate the finite-volume Euclidean matrix element 
{Tnr{n*)\Qi\K)v with the absolute value of the phys¬ 
ical amplitude {'K'K\Qi\K)'. 


{'K'K\Qi\K) = '/t {'K'K{n*)\Qi\K)v , 

^2 pv{E)EmK 


T = 327r^E^ 


AE) 


where 

<E) 

Pv{E) 



dn _ q4>'{q) + k6'{k) 
dE 47rfc^ 


( 10 ) 


( 11 ) 


The last expression in eq. ( 0 ) is the one which one 
would heuristically derive from a naive interpreta¬ 
tion of pv{E) as the density of states, cf. eq. (^). 
There are however, some technical subtleties with 
such an interpretation which will not be discussed 
here. Eq. ( 0 ) holds also when rriK A and the 
operator carries non-zero energy-momentum []l5| . By 
varying the kaon and pion masses and momenta, one 
may then ht the coefficients of the chiral expansion of 
( 7 r 7 r|Qi|itr) and use these coefficients to extrapolate 
to the physical point. This procedure is necessary, at 
present, since it is not possible to compute the matrix 
elements at realistic values of the quark masses and 
in the unquenched case. To give an explicit exam¬ 
ple, the matrix elements of Qr,8 (for generic meson 
masses, one pion at rest and the other with a given 
momentum corresponding to an energy can be 
written as [f^,p5[ 


E^) = 7 ’^’® 

I ^ tti-j^E-jy j - t - ... . ( 12 ) 


A fit to the lattice data for extracted from 

suitable correlation functions, allows the determina¬ 
tion of the couplings 7 ^’®, etc. The results of 

^ For n = 0, there are two solutions: one corresponding to 
k = 0 which is spurious, the other giving the Liischer relation 
between the finite volume energy and the scattering length. 



Figure 1. a. Chiral behaviour of the matrix element of 
< 54 , Ett). The extrapolation to the physical 

point using only the data corresponding to the lightest 
masses, without the effect of the chiral logarithms and us¬ 
ing the operators of ref. [Q, is also shown as a shadowed 
band. The experimental result, computed as explained 
in the text, is also given as a line. b. Chiral behaviour 
of the matrix element of Qg, Mg{mK,Tn.n,E-,r)- The ex¬ 
trapolation to the physical point using only the data cor¬ 
responding to the lightest masses, without the effect of 
the chiral logarithms and using the operators of ref. 
is also shown as a shadowed band. 


the extrapolated amplitudes, if chiral logarithms are 
included, are independent of the cutoff used in the 
chiral theory. 

4. Numerical Results 

In this section, the latest numerical results from 
lattice QCD for K ^ tttt and K ^ tt amplitudes 
and for Bk are reviewed, together with a comparison 
with other approaches. 

A/ = 3/2 transitions and K^—K^ mixing 

The SPQcdR Collaboration has presented 
(quenched) results for both A/ = 3/2 and A/ = 1/2 
amplitudes, obtained from a calculation oi K ^ tttt 
matrix elements following the strategy of refs, [p^ , 
lil- For A/ = 3/2 transitions, an extended study of 
the chiral behaviour of the matrix elements of the op¬ 
erators < 54 , 7,81 renormalized non-perturbatively, was 
performed. In figs. and 0 b, AI 4 (wif, and 

M-sijnKi are shown as a function of the kaon 

mass. Note that the amplitude does depend on three 
independent quantities that are let to vary, namely 
mK^ rriT^ and Et^ (one of the two pions is always at 
rest). In the same figures the extrapolation to the 
physical point, performed by using the chiral expan¬ 
sion of the matrix elements as in eq. (^^, is given as a 
band. In fig. ^ the experimental number, extracted 
from the 7 r+ 7 r° decay rate using the Wilson 

coefficient of <54 computed at the NLO, is also shown. 
The extrapolations are preliminary, since they do not 
include the effects of (quenched) chiral logarithms 
that have not been computed yet for the kinemat- 
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ical configurations used in this study. For some of 
the points, masses and momenta are too large to use 
chiral perturbation theory, and they have not been 
included in the fit. The preliminary results in fig. 
already give us an interesting physics information. 
In the chiral limit, and using SU{3) symmetry, one 
may relate the K~^ amplitude to the K^- 

K° mixing parameter Bk j^ ]- In this limit, a large 
value for Bk, as found in lattice calculations and 
unitarity triangle analyses [^, Bk ~ 0.85, would 
lead to a value of the K~^ amplitude larger 

than the experimental one by ~ 50%. The extrapo¬ 
lation of the results of fig. |^, together with those for 
Bk obtained by the same collaboration given in 
table 0 demonstrate that chiral corrections to both 
the > 7r+7r° amplitude and Bk can easily rec¬ 

onciliate a large value of Bk with the experimental 
^ 7r+7r° amplitude. 

It is very instructive to compare the results for Q^ 
and Qg with those obtained in other approaches, ta¬ 
ble In the case of Qg we notice the nice agreement 
between lattice results obtained from K ^ tt matrix 
elements using chiral perturbation theory and those 
from the first direct calculation oi K tttt ampli¬ 
tudes. Similar numbers were obtained by the CP- 
PACS dl and RBC [|| collaborations, which com¬ 
puted K ^ TV matrix elements with DWF. These 
groups presented the results in different renormal¬ 
ization schemes and scales and for this reason the 
values have not been included in the table. There 
is, however, a very large discrepancy of the lattice 
results with those obtained with dispersive meth¬ 
ods [Q or with the l/N^ expansion whose re¬ 
sults would correspond to a huge value of the B- 
parameter, Rg = 3 ^ 7 |^. In order to reproduce 
the experimental results for e'/e within the Standard 
Model, such a large value of Bg implies a stratosferic 
value for Bq. After Kaon 1999 a new analysis, per¬ 
formed with spectral function techniques by Bijnens, 
Gamiz and Trades appeared In this paper, a 

value of (Qs) much lower than in refs. was 

found, although still about a factor of two larger than 
lattice determinations. The very low value of (Q?) 
found from the lattice K —>■ tttt calculation originates 
from large cancellations occuring in the renormal¬ 
ization of the relevant operator and requires further 
investigation. 

Lattice predictions for Bk have been very sta¬ 
ble over the years, with a central value centered 
around 0.85. This is essentially also the value 
extracted from Unitarity Triangle Analyses [^, 


j_^ther results for Q^ and Qg can be found in refs. [ pH 


and 



Figure 2. First signal for the matrix element (7r7r|Q“|A) 
at the matching point rriK ~ Wo, obtained with 400 con¬ 
figurations and V = 24® x 64. 


thus confirming that lattice QCD can predict (and 
not only postdict) physical quantities. In table 
the lattice world average is given, together with 
some of the latest lattice results. Bk from CP- 
PACS |]^ (with operators renormalized perturba- 
tively) and RBC (with operators renormalized 
non-perturbatively) has been obtained with DWF. 
This formulation of QCD should garantee a better 
control of the chiral behaviour of the regularized the¬ 
ory with respect to Wilson-like fermions. The re¬ 
sults from ref. 1^ have been obtained using the non- 
perturbatively improved Wilson-like action, with a 
new method based on the Ward Identities , which 
save us from the painful subtractions of the wrong 
chirality operators, which was necessary in the past. 
Similar strategies have been pursued with twisted- 
mass fermions [Q. The results with DWF are about 
15% below the world average and show a marked de¬ 
crease at small quark masses. This could reconciliate 
the large value of Bk at the physical kaon mass with 
the low value of this parameter obtained in the chi¬ 
ral limit by ref. [Q. Whether the decrease at small 
quark masses is a physical effect, or is due to lat¬ 
tice artefacts (finite volume, residual chiral symme¬ 
try breaking etc.) remains to be investigated. New 
calculations of {K^\Qi\K°) for all possible operators 
Qi which can mediate K^-K^ mixing in extensions 
of the Standard Model have be presented . These 
results are very useful to put severe constraints on 
FCNC parameters of SUSY models |^|. 

AI = 1/2 transitions and t' je. 

For A/ =1/2 transitions and e'/e, the subtrac¬ 
tions of the power divergencies, necessary to ob¬ 
tain finite matrix elements, are the major obsta¬ 
cle in lattice calculations. These divergencies arise 
from penguin contractions of the relevant operators, 
which induce a mixing with operators of lower di¬ 
mensions, namely sd 

and sysd Q. Power divergencies and subtractions 
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Table 3 

K —> TTTT matrix elements in GeV at ^ = 2 GeV in the MS scheme, from lattice calculations (hrst three rows) and 
other approaches. 


Quantity 


Experiments 

Theory 


Reference 


Method 

(Qs) 

(Qy) 

SPQccR [! 

1| 2001 

K TTTT 

0.53 ±0.06 

0.02 ±0.01 

SPQcbR |E 

) 2001 

K ^ TT + yPT 

0.49 ±0.06 

0.10 ±0.03 

APE PI 1999 

K ^ IT + yPT 

0.50 ±0.10 

0.11 ±0.04 

Amherst |60| 1999 

Dispersion relations 

2.22 ±0.67 

0.16 ±0.10 

Marseille p 

l| 1998 

1/N, 

3.50 ± 1.10 

0.11 ±0.03 

BGP [|^ 2001 

Spectral functions 

1.2 ±0.5 

0.26 ±0.03 


Table 4 _ 

Bk in the MS scheme at the renormalization scale /r = 2 GeV and the renormalization group invariant Bx 
from recent lattice calculations. The results from SPQcdR., CP-PACS and RBC Collaborations are quenched. 


Quantity 

Experiments 

Theory 


Reference 

Method 

B^^{2 GeV) 

Bk 

World Average 

SPQcuR H 2001 
SPQcnR ^001 
GP-PACS g 2001 

RBG [H 2000 

Giuchini et al. 2001 

with subtractions 

Ward identity method 
DWF Pert. Ren. 

DWF Nonpert. Ren. 
Unitarity Triangle 

0.63 ±0.03 ±0.10 
0.71 ±0.13 

0.70 ±0.10 
0.5746(61)(191) 
0.538 ±0.08 

0 70+®-^® 

0.86 ±0.06 ±14 
0.91 ±0.17 

0.90 ±0.13 

0.787 ±0.008 
0.737 ±0.011 

0 nn+0.30 


6.Qa-oi I--— 

4.0S-01 - 

2.0S-01 - 

o.os+oo ' ‘ • I j j J 
-2.0S-01 - 


Figure 3. First signal for the matrix element {mT\Qe\K) 

at the matching point rriK ~ Wd, obtained with 400 con- 

hgurations and V — 24® x 64. 

6.06-01 

4.06-01 

2.06-01 

I 

0.06+00 

-2.06-01 


Figure 4. (7r7r|s75(i|A') at the matching point rriK ~ 

Wo, obtained with 400 configurations and V = 24® x 64. 
At the matching point (7r7r|s75d|A) decreases by a huge 
factor and thus the subtraction for {-kttIQoIK) is very 
small. 


are encountered with both the strategies (AT ^ tttt 
or A ^ tt) adopted to compute the physical am¬ 
plitudes. For example, with a propagating charm 
quark. Re Aq is computed in terms of the matrix el¬ 
ements of {Qi and Q 2 ) only. In this case, due 
to the GIM mechanism, the subtraction is implicit 
in the difference of penguin diagrams with a charm 
and an up quark propagating in the loop: 

= S7^(l-75)'uu7'"(l-75)d (13) 

± 57^ (1 - 75) dM7'^ (1 - 75) u - (m ^ c). 

In the past any attempt to compute AI =1/2 
K —>■ TTTT matrix elements failed because no visi¬ 
ble signal was observed after the subtraction of the 
power divergencies These calculations, per¬ 

formed about twelve years ago with very modest 
computer resources, on small lattices and with lit¬ 
tle statistics, were abandoned after the publication 
of the Maiani-Testa no-go theorem . 

This year, the SPQcuR. Gollaboration has found 
the first signals for both the matrix elements of Q~ 
and Qq, figs. |, | and i The difference with re¬ 
spect to previous attempts is given by a much larger 
volume and statistics, about 400 configurations with 
V = 24® X 64, compared to a few tens of configura¬ 
tions on lattices with V = 16® x 32 (8 configurations 
on V = 24® X 40) or 110 configuration of a lattice 


t i n H H f n H H i J J j 










































with V = 16 X 12 X 10 X 32 (sic !!) (Q. A crucial 
ingredient is also to work at the point correspond¬ 
ing to the matching condition mx = Wq, where 
Wq is the energy of the two pions at rest on the fi¬ 
nite volume used in the simulation, namely with non 
degenerate strange and light quark masses Past 
calculations were always performed, instead, at the 
degenerate point, rriK = rriTj. The results are very 
preliminary and in particular the two-pion energy in 
the finite volume for this channel gives a scattering 
lenght, a/=o, in disagreement with expectations both 
in value and in the mass dependence. This is to be 
contrasted with the A/ = 3/2 case, where the analy¬ 
sis of the scattering length ai =2 is in good agreement 
with expectations |^. Much more work is needed to 
clarify these problems before trying to extrapolate 
the amplitude to the physical point. 

At the matching point niK = Ho, there are also 
results for {■mT\QQ\K) . In this case there is no GIM 
mechanism at work, and a finite subtraction of the 
matrix element of the pseudoscalar density opera¬ 
tor Qp = (rus — md)s')^d must be done (for par¬ 
ity, instead, the subtraction of the scalar operator 
Qs = (ws -|- md)sd is not necessary). The coeffi¬ 
cient of the mixing of Qp with Qg is quadratically 
divergent ^ 

4 Qe + (14) 

a‘‘ 

If not for the explicit chiral symmetry breaking of 
the lattice action, {'ktt\Qp\K) would vanish by the 
equation of motions when rriK = Wn* ■ With an 
improved action, {-k'k\Qp\K) is of 0{a^) and thus 
a finite subtraction is necessary. The preliminary 
results shown in fig. ^ show that that the value of 
{'KTT\s^^d\K) drops by a factor of about 20 with re¬ 
spect to {'K'k\Qq\K) when ttik = Wq (it is of the same 
order when the kinematics is not matched) and thus 
the subtraction is not very critical. The signal itself 
is very noisy however, see fig. |^, and even with 400 
configurations the statistical error on {'ktt\Qq\K) is 
about 50%. Thus a rather large statistics will be 
necessary to obtain a relatively accurate result. 

K TTTT amplitudes in the chiral limit can also 
be extracted from the calculation of AT ^ tt ma¬ 
trix elements, and indeed this has been the most 
popular method in lattice calculations of A/ = 3/2 
transitions. The main advantage of this approach is 

® In this case, as in the case of degenerate quark masses, 
the subtraction of operators of lower dimensions is not nec¬ 
essary ii- 

The mixing with the chromomagnetic operator is not dis¬ 
cussed here. This mixing is a small, finite correction which 
can be ignored to simplify the discussion. 


that the three-point correlators, necessary for the ex¬ 
traction of the matrix elements, are less noisy than 
the four-point correlators used in the K tttt case. 
Moreover, for K ^ t: matrix elements, finite volume 
corrections are exponentially small as L oo. The 
main disadvantage is that the K tttt amplitudes 
are obtained, using soft pion theorems, in the chiral 
limit only and then the effect of FSI is definitively 
lost. If these are important for A/ =1/2 channels 
there is no hope, then, to recover the physical ampli¬ 
tudes. In order to make the IF —> tt matrix element 
finite at least a subtraction is necessary, even in the 
absence of explicit chiral symmetry breaking in the 
action. With Wilson-like Fermions, for which chiral 
symmetry is explicitly broken, the number of sub¬ 
tractions makes this approach unpracticable as also 
demonstrated by the failure of past attempts @11 • 
The method has acquired a renewed popularity, in¬ 
stead, with recent formulations of the lattice theory, 
DWF or overlap Fermions, for which chiral symme¬ 
try breaking is absent or strongly reduced in practice. 
In this case, and under the hypothesis that chiral 
symmetry is under control, a finite amplitude can be 
obtained by subtracting the scalar density amplitude 
with a suitable coefficient Ci which depends on the 
operator at hand 

{nlQr^lK) = {nmK) - C,{7r\Qs\K). (15) 

The power divergent coefficients Ci cannot be com¬ 
puted perturbatively. Following ref. , when chiral 
symmetry is preserved, the Ci can be determined by 
the condition that the K ^ 0 matrix element of the 
subtracted operator vanishes 

{0\Q^-C,Qp\K) = 0. (16) 

The coefficients Ci have been obtained either using 
non degenerate quarks, rUs ^ irid, by the RBC Col¬ 
laboration, or from the derivatives of the 2-point cor¬ 
relation function with respect to the quark mass, 
by the CP-PACS Collaboration. Numerically, the 
matrix element of the subtracted operator is much 
smaller than the unsubtracted one, as shown by 
fig-H taken from the recent work of the RBC collab¬ 
oration [Q. This implies that any systematic uncer¬ 
tainty which enters in the subtraction procedure can 
have huge effects in the determination of the physical 
amplitudes. After the subtraction, the chiral depen¬ 
dence of the matrix element has to be fitted in or¬ 
der to extrapolate it to the chiral limit. Both groups 
have included the logarithmic corrections which arise 
in quenched xPT in the fit and, in some cases, poly¬ 
nomial corrections of higher order in mj. The chiral 
behaviour observed by CP-PACS is very satisfactory, 
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-0.01 0 0.01 0.02 0.03 0.04 0.05 


Figure 5. {n\Q\K) for the unsubtracted operator, Qq 
(squares), the subtraction, C&Qs (Circles), and the total, 
Qf“*’ see eq. di), as a function of the quark mass m/. 
The results are from the RBC Collaboration. 



[GeV’l 


Figure 6. Chiral behaviour of the matrix element 
{'k\Qq\K) from the CP-PACS Collaboration. The curves 
represent different fits used to extrapolate to the chiral 
limit. 

as shown by fig. less good in the case of the RBC 
Collaboration. The difference may be due to the fact 
that a different gluon action is used in the two cases, 
corresponding to smaller chirally breaking effects for 
CP-PACS. This can be monitored by measuring the 
so called “residual mass” which should be zero for 
perfect chirality and is about a factor of ten smaller 
for CP-PCAS than for RBC. For a more extended 
discussion on this important point the reader can re¬ 
fer to Let me also mention that the overall 

renormalization constants have been computed per- 
turbatively by CP-PACS Q and non-perturbatively 
by RBC [^. My compendium of the physics results 
obtained by the two groups is given in table A 
few observations are necessary at this point. First of 
all, i?eMo, and consequently ReAo/ReA 2 , from the 
RBC Collaboration are in good agreement with ex¬ 
perimental values, contrary to CP-PACS which finds 
RcAq smaller by about a factor of two than the ex¬ 
perimental number. Since the two groups use the 
same lattice formulation of the theory and differ only 
by, hopefully, marginal details, the reason of this dif¬ 
ference should be clarified. In both cases, however, 
e'/e is in total disagreement with the data. The main 
reason is that the value of the matrix element of Qq is 
much smaller than what would be necessary to repro¬ 


duce the experimental value (it corresponds approx- 
imatively to Bq = 0.3-b0.4). Let me list a number of 
sources of systematic errors which may explain these 
embarassing results: 

1) Chirality By working with DWF at a finite 
fifth dimension, TVs = 16, a residual chiral symme¬ 
try breaking is present in the theory. The amount 
of residual symmetry breaking is parametrized by 
a mass scale denoted by mres 0.2 -h 2.0 MeV. 
In the presence of explicit chiral symmetry break¬ 
ing, the coefficients Ci determined from K ^ 0 dif¬ 
fers from the correct one and this may induce an 
error of 0{mresa~^) ~ (200 MeV)^ on matrix ele¬ 
ments which are of the order of ^qcd ~ (3*^^ MeV)^. 
Both groups claim to have this point fully under con¬ 
trol H,^. A calculation at a larger value of N 5 , with 
all the other parameters unchanged, would be very 
useful to clarify the situation. 

2) Matching below the charm mass Both groups 
have so far presented results at a renormalization 
scale just below the charm quark mass. The match¬ 
ing of the effective theory is rather problematic at 
such low scales. 

3) Extra Quenched Chiral Logarithms As shown 
by figTir the subtraction of the power divergencies 
is very critical. Besides the effects discussed in 1), 
another delicate problem has been recently rised by 
Golterman and Pallante 113|. In the quenched theory 
the (8,1) operators, such as Q~ and Qe, do not be¬ 
long anymore to irreducible representations of the 
chiral group and this gives rise to spurious chiral 
logarithms. These logarithms affect the subtraction 
procedure and have not been taken into account in 
the analyses of RBC and CP-PACS. A preliminary 
numerical study to try to get rid of these effect has 
been discussed at this conference 

4) Higher order chiral corrections and FSI Higher 
order chiral contributions, among which FSI have 
also to be accounted, can produce large variations 
of the matrix elements between the chiral limit and 
the physical point. In simulations where only K tt 
matrix elements are computed, K —> tttt amplitudes 
can only be obtained at lowest order in yPT where 
these physical effects are missing. 

5) New Physics As noticed by Murayama and 
Masiero nr it is possible to produce large effects for 
e'/e in SUSY without violating other bounds coming 
from FCNC. If the lattice results are correct, this 
is an open possibility. Before invoking new physics, 
at least the result for ReAo should be established 
with more confidence. Given that the two groups still 
don’t agree on this quantity, more work is needed. 
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5. Conclusion and outlook 

In this review, because of space limitations, I dis¬ 
cussed only a subset of the results for light-quark ma¬ 
trix elements presented at Lattice 2001. Although 
the focus has been on kaon weak matrix elements, 
the general discussion relative to renormalization, 
improvement and calculation of physical amplitudes 
applies to many other cases as well. A renewed activ¬ 
ity in lattice calculations of kaon decays and mixing 
has developed in the last two years. For AI = 3/2 
transitions, K'^ —> and electroweak penguin 

amplitudes, new and more precise results have been 
obtained. In this case, removal of lattice artefacts, 
by extrapolation to the continuum and/or improve¬ 
ment, and unquenched calculations are around the 
corner. For AI = 1/2 transitions, it has been shown 
that direct computations of K tttt amplitudes, 
including FSI, is possible (0111, although the prat- 
ical implementation with sufficient accuracy will re¬ 
quire more time. First results, obtained by using 
AT —> TT lattice matrix elements computed with DWF 
and soft-pion theorems, show a striking disagreement 
with the experiments for e'/e. For RcAq, there is a 
factor of two between the values found by the CP- 
PACS Q and RBC ||] Collaborations. More work is 
needed in this case to clarify all these points. Many 
related physical quantities, like the strong interaction 
phase-shifts for tttt scattering, the chromomagnetic 
operator matrix elements, semileptonic K —>■ nnli'i 
amplitudes and the scalar semileptonic form factor 
will also be obtained as a byproduct from lattice 
investigations of non-leptonic kaon decays. In par¬ 
ticular, the study of the K ■k'kIui amplitude is 
very useful since it contains all the ingredients rela¬ 
tive to FSI without having the extra difficulties due 
to the power divergences of the non-leptonic ampli¬ 
tudes, and the corresponding necessary subtractions. 
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